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Abstract 

In this paper, we show that an attempt to construct shape invariant extensions of a 
known shape invariant potential leads to, apart from a shift by a constant, the well known 
technique of isospectral shift deformation. Using this, we construct infinite sets of gen¬ 
eralized potentials with Xm exceptional polynomials as solutions. These potentials are 
rational extensions of the existing shape invariant potentials. The method is elucidated 
using the radial oscillator and the trigonometric Poschl-Teller potentials. For the case of 
radial oscillator, in addition to the known rational extensions, we construct two infinite 
sets of rational extensions, which seem to be less studied. For one of the potential, we 
show that its solutions involve a third type of exceptional Laguerre polynomials. Explicit 
expressions of this generalized infinite set of potentials and the corresponding solutions 
are presented. For the trigonometric Poschl-Teller potential, our analysis points to the 
possibility of several rational extensions beyond those known in literature. 
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1 Introduction 


The search for new exactly solvable (ES) models gained momentum with the discovery of 
the exceptional orthogonal polynomials (EOPs) [I],[l]- New potentials with the Xi EOPs as 
solutions were constructed in mM- They turned out to be the rational extensions of the well 
studied one dimensional bound state problems namely the radial oscillator and the Darboux 
Pdschl-Teller (DPT) potentials [5]. A spate of papers followed, where generalized rational 
extensions of many one dimensional shape invariant potentials (SIPs) were presented [5]- |17j . 
Different methods like the Darboux-Crum method finite difference Backliind algorithm 

m-m, prepotential method m and several others m, m were employed. Various groups 
explored the mathematical properties of the new polynomials concurrently m-m- 

Many studies showed that the traditional SIPs and their rational extensions are isospectral 
and share a supersymmetric partnership 0, m, m In recent papers, Grandati et.al., |12j . 
m constructed rational extensions of SIPs using regularized excited states Riccati-Schrodinger 
functions as superpotentials. The finite difference Backlund transformation was combined with 
the regularization scheme to generate infinite sets of rational extensions of the isotonic oscillator 
|18) . Dutta, Roy and coworkers [16], [21] and references therein, constructed new conditionally 
exactly solvable (CES) models using supersymmetric techniques. Among these a particular 
construction gave potentials related to exceptional polynomials Ell- 

Shape invariance (SI) property, discovered in [^ , has been central to the discussion of the 
exact solutions of potential problems in supersymmetric quantum mechanics (SUSYQM) j^, 
m, m- We make this property the focus of our study in this paper. Here we formulate a 
sufficient condition for SI using the quantum Hamilton Jacobi (QHJ) formalism - [30]. This 
condition brings out the fact that SI is a property of a potential and does not require knowledge 
or a reference to a SUSY partner potential or the superpotential. This requirement makes it 
possible to analyze the SI property and use it to obtain deformed potentials which are shape 
invariant by construction. This will be illustrated using the radial oscillator and the Darboux 
Pdschl-Teller (DPT) potentials. In an earlier paper [31], it has been shown that all previously 
known SIPs, as listed in [^, can be derived in an extremely simple, straightforward method using 
the sufficient condition for SI. In this paper we complete this program by constructing rational 
SIPs of radial oscillator leading to potentials related to the LI and L2 exceptional polynomials 
0, m. [H], m, M- In addition, we also find one more generalized set of potentials with 
polynomial solutions, which may be classified as type HI or L3 exceptional polynomials m,m- 
The expressions for these generalized potentials and their polynomial solutions are presented 
here for the first time. Eor the DPT potential, we are led to many rational extensions of which 
some coincide with the known extensions in the literature 0 , 0 , m- 

Our scheme allows us to construct a hierarchy of generalized ES rational potentials indexed 
by the hierarchy index m, (to = 0,1, 2...). Each value of to gives a new rational extension of 
the original SIP and the solutions are in terms of the EOPs. For to = 0, the generalized 
potentials reduce to the original SIPs. In the scheme presented in this paper, the superpotentials 
of the original SIP, corresponding to both exact and broken SUSY, play a central role. The 
interplay of ideas and the use of tools available in the QHJ formalism and SUSYQM makes our 
scheme very simple, elegant and transparent. 

The knowledge of the singularity structure of the quantum momentum function (QMF) 
within the QHJ formalism, allows one to construct all possible superpotentials associated with 
a given supersymmetric potential. Next for each superpotential, say W{x)^ we find extended 
superpotential, W{x), such that the corresponding potentials U^“''^(a;) is shifted w.r.t. U(+)(x) 
by a constant. Each of the resulting potentials, U(“^(a;), associated with W{x), is either a 
singular or a regular rational extension of V^~\x). This method has been earlier used in 
m to obtain CES rational potentials with EOPs as solutions. The superpotentials associated 
with nonnormalizable solutions of V^~\x) lead to the regular extensions and these potentials 
have eigenfunctions in terms of the Xm exceptional polynomials as solutions. The potentials 
associated with the normalizable wave functions result in new sets of generalized potentials 
which are singular or regular depending on the potential parameters. Some of these potentials 
are discussed in 0, [10 . 

In Sec. 2 we first recall some important features of the QHJ formalism of quantum mechanics 
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and reformulate the SI requirement, crucial to our discussion. In Sec.3, we show that this 
requirement leads to the isospectral shift deformation of the SIPs known from the early days of 
SUSYQM. Sec.4 has a detailed discussion of rational extensions of the radial oscillator. We also 
present the expressions for an infinite set of rational potentials, whose solutions are in terms of 
the new set of polynomials, the L 3 category of exceptional Laguerre polynomials. The case of 
trigonometric DPT, leading to the exceptional Jacobi polynomials is taken up in Sec.5. Finally 
we present our conclusions in Sec. 6 . 


2 The QHJ and SUSYQM formalisms 

In the QHJ formalism - [^, the starting point is the Riccati equation 

E) — ih ’ —^-= 2m{E — V(x)), (1) 

ax 

obtained by making the substitution ip{x,E) = exp(iS{x, E)/H) in the Schrodinger equation. 
Here p{x,E), the quantum momentum function (QMF), is identified with ■ For conve¬ 

nience, we will use units so that h = 1, 2m = 1 and in order to work with real equations, we 
define Q{x, E) by p{x, E) = iQ{x, E). The QHJ equation for Q{x, E) assumes the form 


Q\x,E) 


dQ{x, E) 
dx 


V{x) - E. 


( 2 ) 


Q{x, E) is simply the negative of the logarithmic derivative of solution E) of the Schrodinger 
equation 


Q{x,E) 


jL 

dx 


\ogtlj{x,E) 


1 d'tjj{x,E) 
'0(x, E) dx 


(3) 


The QHJ approach essentially makes use of the singularities of the function Q{x, E) in the 
complex X plane. A crucial step in solving for the QMF, Q{x, E), is to write it as a sum of two 
parts, 

Q(x, E)=Qf + (4) 


where Qf contains the fixed singularities and its form can be written down by inspection of 
different terms in the QHJ equation. The second part, Qm, contains the singular part having 
the moving singularities which depend on the boundary conditions to be imposed. In all known 
cases of ES potentials Ez], [28] including the new rational potentials with EOP solutions [29] and 
the quasi-exactly solvable potentials [50]) when E is taken to be one of the energy eigenvalues, 
Qm turns out to have a finite number of simple poles with residue -1. It is easy to see from m 
that the moving poles of Q{x) are located at the nodes and at other complex zeros of the wave 
function. 

The superpotential w{x), which in SUSYQM corresponds to the negative of logarithmic 
derivative of a solution of the Schrodinger equation, plays an important role of defining the 
partner potential and formulating the SI condition [52]. Usually, the solution chosen has no 
nodes but may, or may not, be normalizable, the two cases being referred to as cases of exact 
and broken SUSY respectively. Thus w{x) happens to be just one of the solutions Q{x, E) 
of the QHJ equation, for certain energy, for which moving poles are absent and Qm in 0 is 
replaced by a constant. The SUSY partner potentials V^{x) are defined by 


= w{x)'^ ± w'{x) 


(5) 


and a potential V(x) is defined to be shape invariant if both the partner potentials have the 
same form as V(x) apart from an overall constant. Since several superpotentials are possible, 
hence the definition of partners, and the property of SI itself appears to depend on the choice 
of the superpotential. It turns out to be extremely useful to formulate SI as a property which 
does not require the use of the superpotential, and which makes it transparent as a property 
of the potential. As we will see later, this fact assumes importance in the context of rational 
extensions related to the exceptional polynomials. 
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We now reformulate the SI requirement on a potential V(x). Let Q{x,a) be a solution of 
the QHJ equation, 

Q^(x,g) - =V{x) - E{a), (6) 

where the dependence of Q{x, a) on the potential parameters has been explicitly shown through 
constants cr, which are functions of the parameters appearing in the potential. If we introduce 
two functions wi(x), W 2 (x) 

wi(x) = Q(x, A), W 2 (x) = Q(x,fi). (7) 

corresponding to two different values cr = A, then the two potentials, defined by 

Vi(x) =Wi(x) -w[(x), V 2 (x) = wi(x) - W 2 (x), (8) 

obviously will have the same shape. If we further require that there exists a map 

T : X ^ ^ = t(A) (9) 

such that 

Q{x,X) +Q{x,n) =0, (10) 

then W 2 {x) = —'Wi{x) and the two potentials Vi, 2 (a^) will be SUSY partners. 

To summarize, given a potential V{x), the twin requirements of existence of solutions 
Q{x, A), Q{x, pf) of QHJE and of a map r such that 

Q{x,n) =-Q{x,t{X)) (11) 

are seen to be the sufficient conditions for SI of the potential V (x). 


3 Shape Invariant Rational Extensions 

Given a superpotential w{x, A), the SI requirement means the potential V^'^\x, A) = w‘^{x, A) + 
w'{x,X) be equal to V^~\x,p) = w‘^{x,p) — w'{x,p), apart from an additive constant, after 
the set of potential parameters. A, is redefined. This translates into the following requirement 
on the superpotential as 

w‘^{x, X) + w'{x, A) = w'^{x, p) — w'{x, p) + constant. (12) 

for some A, p and a constant which may depend on A, p. 

Let us now assume that a particular superpotential wq{x, X) has been found which gives rise 
to a SIP V{x,X), where A denotes constants which are functions of the potential parameters. 
We now seek to construct a shape invariant extension of wo(x, A) by constructing an extended 
superpotential 

w(x, X) = wo(x, X) + f)(x, X), (13) 

where (/)(x, A) is an unknown function and is to be determined from the requirements of SI. 

Using the fact that both w(x) and wo(x) for some p, A, satisfy the SI requirement (HU , we 
get 


(j){x, X)^ + 2wo{x, X)4i{x, A) + (f'ix, A) (14) 

= (j){x,p)^+ 2wo{x,p)<p{x,p) - (j)’{x,p) + K, (15) 

where K is a, constant. For SIPs of interest in this paper, and defined by wo{x, A), there exists 
a map t : X ^ p = 'r(A) such that 

woix, p) = -wo{x,t{X)). (16) 

We also want the extended potential to be shape invariant and hence we demand that 

(l){x,X) ^ -(I){x,t{X)). (17) 
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Thus we can rewrite (fT^ as 


(j){x, X)^ + 2wo{x, X)(j){x, A) + 4‘'{x, A) (18) 

= (j}{x, r(A))^ + 2wo{x, t{X))4>{x, t(A)) + (j)'{x, r(A)) + K. (19) 


The structure of the above equation suggests using the ansatz of equating both sides of the above 
equation to a constant. Therefore, we look for solutions for 4’{x, A) satisfying the equation 

(/)(x, A)^ + 2wo{x, X)4>{x, A) + (j)'{x, A) = Ri, (20) 


where Ri is a constant depending on the potential parameters A. The above condition, (1201), is 
simply the relation 

V'^+\x) = V^+\x)+Ri, ( 21 ) 

where the notation y(=*^)(a;) and lA^^)(x) denote the SUSY partner potentials constructed from 
w{x) and wo(x) respectively. It may be remarked that we have recovered the property that has 
been used in m to construct new conditionally exactly solvable models with solutions related 
to the exceptional polynomials. For the special choice = 0, this requirement is same as that 
used in construction of strictly isospectral deformations of a given potential [32] , [23] , 

U(“)(a:) = w;o(a:, A) - u;o(a;. A). (22) 


In general the constant Ri in 
by means of 


(EU can be non zero and the above process of defining ) (x) 


^x) = w'^{x,X) — w'{x,X) 


(23) 


will be called the isospectral shift deformation of potential U^“)(a;). 

The equation for (j){x, A) is the Riccati equation and can be linearised by introducing u(x, A) 
defined by 


(j){x, A) 


1 du{x, A) 
u{x, A) dx 


(24) 


The unknown function u{x, A) then satisfies a Schrodinger like equation 


d^u{x, A) 
dx'^ 


+ 2wo{x, X)u{x, A) — Riu(x, X) = 0. 


(25) 


In order to obtain rational extensions and to connect with the exceptional polynomials, the 
constant Ri will be chosen so as to give polynomial solutions for u{x). Let Pm{x) denote a 
polynomial solution, of degree m, of (1351) . A rational extension of the original SIP V{x), will 
be obtained from the polynomial solutions Pm[x) for u(x) by means of the relations 


^(x) = w^(x) — w'(x), 


(26) 


where 

1 f] P ( 

w{x)=wo{x) + cf{x), cl>{x)= (27) 

The equations (1331) - (1371) are the basic equations for the rational extensions and to arrive at 
the Hamiltonians which have eigenfunctions in terms of the EOPs. 

In general when Ri is not restricted and is allowed real values, we would get an extended 
potential which would interpolate between the potentials related to the exceptional polynomials. 


4 Rational Extensions of the Radial Oscillator 


4.1 Solutions for the Superpotential 

In this section we will construct the rational extensions of the radial oscillator given by ( 2m = 1) 


V{r) = + 




(28) 
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Four superpotentials w{r) for the radial oscillator, obtained from the solutions of the QHJ 
equation, 

w^{r) — w'{r) = V{r) — E, (29) 

and the corresponding partner potentials (r) are given in the table below. 


Table 1 


k 

Superpotential Wk 

II 

s 

1 

yt'' =wl + w'k 

1 

1 {£+!) 

-Lur - 

2 r 

V{r) -ujie + 3/2) 


2 

1 i 

-ur + - 

2 r 

V(r)+uji£-l/2) 


3 

1 (£ + 1) 

—-cjr- 

2 r 

V{r)+uj{£ + 3/2) 


4 

1 i 

--wr + - 

2 r 

V{r) -uj{£-l/2) 

\uV + £l^-u,(t+ll2) 


These four solutions of the QHJ equation are meromorphic, have no moving poles and are 
of the form ^ 

w{r, cr) = - + -br, cr = {a, b}, (31) 

where a, b have two possible values 

a = -(£ + !),£; b = ±uj. (32) 

The four superpotentials correspond to the four possible combinations of values for {a, 6}. The 
parameter values, 

fi = {-£-l,u;}, X = {£,-u)}. (33) 

correspond to the two solutions < 1 > and < 4 > in (1501) of QHJE. 

With translation mapping r defined as 

'r({^)W}) = {£+ l,w} 

we have t(A) = {£+!, —uj}. The SI property (ITTl) . 

wir,fJ-) =-wir,T{X)), (34) 

for the radial oscillator, is obviously satisfied. The solutions W 2 {x), W 3 (x) for the superpotential, 
in the second and the third rows of the table correspond to broken supersymmetry. These also 
are related by a relation of the form (IMll . These will be used for constructing rational potentials 
related to the Li and L 2 type exceptional Laguerre polynomials. 

The extended potentials corresponding to wi(r) and 104 (r) are also constructed. Most of 
these potentials turn out to be singular. However in some cases regular rational extensions are 
found and the corresponding polynomial solutions are presented in Sec 4.4. 

4.2 Isospectral Shift Deformation 

The isospectral shift for the radial oscillator will be used to construct a new potential V (r) by 
making use of Eqs. (I211)-(EI1). Identifying wo(r) with W 2 {r), we introduce 

W 2 (r) = W 2 ir) (35) 
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( 36 ) 


and write (j){r) = logu(r), the equation ( 1 ^ then becomes 

/ 2 i\ 

u"{r) + H- ju'{r) — Ru{r) = 0. 

That this equation reduces to the Laguerre equation, is most easily seen by making a change 
of variable r] = — which transforms the above equation to 

r]u"{r]) + (^-r] + i+ = 0 - (37) 

Demanding that the solution for u{r]) be a polynomial in 77 of degree to, gives R = 2muj and 
the polynomial solution coincides with the associated Laguerre polynomials [34], |35| and we 
have 

u( 7 y) = ( 77 ), </)(r) = a = {i- 1/2). (38) 

Therefore, indicating the dependence of the solutions on to, the hierarchy index explicitly, we 
write 

W 2 .m(y) = W 2 (y) ( 39 ) 

where y = = —77 is the variable that will appears in all final expressions. Substituting 

(IMll in (E51) gives the final expression for the rational extension as 


Vi-Hr) 


\r) + 2 uir 


21 

(jjr -\ - 

r 






L^(.y) . 



dyL^^jy) 


2 


(40) 


Idn \r) to show the dependency on the hiararchy index to. For to = 0 all the rational expres¬ 
sions reduce to the undeformed potentials. For each value of to > 0, we get a distinct rational 
extension of the radial oscillator. The partner potential y(+)( 7 -) 


’(/(+) (r) = -|- 2 tow, 

4 


(41) 


by construction, coincides with the radial oscillator. The eigenfunctions of Idn \'r) can be 
obtained by using the (broken) SUSY intertwining relation between the solutions of the partners, 
which for this case becomes 

Yk~L{y) = (^--^+W 2 m{r)^ (42) 

where 

= y^/2exp(-77/2)L“(7/)|j^^i^^2, (43) 

are obtained by replacing £ ^ — 1 in the eigenfunctions of the radial oscillator (110411 in the 

appendix. On simplification, as explained in the appendix, (1421) gives 


= 


exp(—iwr^) 

L^i-y) 


(LL“+'^(-2/)in(2/)-i“(-2/)5,i( 


J y—^uiT'^ 


(44) 


The recurrence relation 

|;L“(r) = T“(r)-T“+i(r) (45) 

has been used to write the eigenfunction in the form (H^ . The details of this step can be found 
in the appendix. The eigenfunctions of the extended potential have the form 




rU 2 


exp ( - 
^mi y^ly—^ujr'^ 


Pn,m ('^); 


(46) 
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where 


Pn,m (^) — 


^\y)-L^ '^\-y)dyL^^ '^\y) 


(47) 


The polynomials Pn,m{'r') are the Xm, EOPs, orthogonal with respect to the weight function 


W(r) = exp 


W 2 ,m(r)d 


exp (— 


(48) 


in the interval 0 < r < oo. The above polynomials have n + m zeros, where n real zeros are 
inside the interval of orthogonality and m zeros, either real or complex, lie outside the interval. 
For all the states, the number of zeros outside the orthogonality interval remains fixed, while 
the number of real zeros located in the interval is governed by the oscillation theorem. These 
polynomials in (|T7| reduce to ’undeformed polynomials’ Ll^^^^(a;r^/ 2 ) when m is set equal to 
zero. From the above, one can see that each value of m gives a different rational extension of 
the radial oscillator (l2^ , whose excited state wave function can be obtained from (1461) . For 
example, putting m = 1 in the equations (j40p and simplifying gives 


V/ \r) = + 


£{£ + 1) + 4a;£ - 2a; 

n 


a;r2 + 2^+1 


1 

*+ 2 £+l )2 +(^- 2 )^^ 


(49) 


which matches with the potential constructed by Quesne in [3]. solutions. Similarly substituting 
m = 1 in (PI). we obtain 




/r^+iexp(-ioar2)\ 

a;r 2 + 2 £ + 1 j 


(ujr^ + 2i+ l)Ll^+^\ia;r2) + 


(50) 


where the polynomial part coincides with the Xi exceptional Laguerre polynomials of Li series 

0 , m- 


4.3 Shape invariance of the extended potential 

We arrived at the new potential V^~\r) by analysing the SI condition (IT^ for the extended 
potential. However, at this stage the SI of the extended potential, though guaranteed, is not 
transparent. This is because the SUSY partner of V^~'^{r) obtained using W 2 ,m is, by definition, 
same as potential Ui“'"^(r) apart from a constant. In order to explicitly demonstrate the SI, 
we obtain a new superpotential satisfying the twin properties given at the end of Sec, 2. For 
this purpose, we now start with W 3 {r) and carry out a similar process with roles of partners 
exchanged. Specifically, we construct a rational extension W{r) = wsir) + x(r) and carry out 
isospectral shift deformation, exchanging roles of U*^^^(r). In this case, the solution for x(r) is 
obtained by demanding 

V^~\r) = V^-\r) + R 2 (51) 

and construct an extended potential 

V^^\r) = w{r)^+w'{r), (52) 

which turns out to be the SUSY partner of V (r) obtained in the previous section. Going 
through the same steps as above, the equation for wir) turns out to be 

X^{r) + 2w3{r)xir) - x'ir) = R 2 , (53) 

where i ?2 is a constant depending on the potential parameters A and 

—/ \ / \ / N / \ 1 dv{r) 

w{r) =W3(r) + x{r), xir) = - 

v(r) dr 

The equation satisfied by v(r) can be seen to Laguerre equation by changing variable to rj = 


Wiv) - -(^ + i) + 


(55) 












Again restricting to polynomial solutions and using y = —rj, we get 


Ky) = 


y=^uir^ 


■j(r) = w^ir) — ujr 


Substituting w{r) in (I52|) gives the expression for (r). We will now find a new superpotential 

Wo{r) such that the two extended potentials V^~\r) and V^^\r) can be seen to be SUSY 
partners, viz., 


a r(«+l) 


i-y) 


L 


(a+1) 


i-y) 


(56) 


U(-)(r) = w2(r)-lU'(r) 
and = lUp (r) + lUo(r). 


(57) 

(58) 


For this we collect all the equations relating the potentials U*^=*=^(r) and V^^\r) and list them 
below. 


U(+)(^) = v'^+'>{r) + Ri, (59) 

F^-\r) = U(-)(r)+i?2, (60) 

U^“^(r) = w{r)'^—w{r), (61) 

l/^'''^(r) = w{r)‘^ +w'{r). (62) 

Note that (15^ implies 

A) + 2w2ir, \)(t)ir, X) = Ri - (j)'{r, A). (63) 

Similarly, the equation for (IHOll implies that x('r) obeys a relation 

X^(r, A) - 2u;3(r, A)x(7’, A) = i?2 + x'(u A). (64) 


Using (l571) in the l.h.s of (ICT|) and substituting tUo(r') = wi{r, A) +'C('r, A) gives w{r) = W 3 {r) + 


(/)(r) to get 

L.H.S. of (ICTI) = i'Wi{r, A) + A))^ — iw[{r, A) + ^'(r, A)) ( 65 ) 

= wlir, A) - w'Fr, A) + ^(r, A)^ + 2 u;i(r, A)f (r, A) - f'(r-. A) 

= U(r) - u;{i + 3 / 2 ) + ^(r, A)^ + 2wi{r, A)'e(r, A) - C'(r, A). ( 66 ) 

Next, the r.h.s of (EH) is expressed in terms of w{r) = W2(r) + (p(r) gives 

R.H.S. of dnH) = (w 2 (r, A) + A) - (^(u A) + </') (67) 

= u'Ku A) - W 2 ir, A) + (/^(r, A) + 2'W2ir, A)(/(r, A) - (/'(r. A) 

= U(r) + uj{£ - 1 / 2 ) + (/^(r, A) + 2w2ir, A)(/(r, A) - A). ( 68 ) 

Equating the expressions in (IM|) and (l 68 )) we get 

C^(U A) + 2 wi(r, A)C(r, A) - ^'(r, A) ( 69 ) 

= 2 a;(£ + 1 ) + (/^(r. A) + 2 w 2 (?', A)(/(r, A) - (/'(r, A) ( 70 ) 

= 2 w(^+l)+i?i- 2 (/'(r-,A), ( 71 ) 

where, in the last step, ( 1551 ) has been used. Similarly from ( 1551 ) and (IMl) we would get 

^ 2 (r, A) + 2 u;i(r, A)C(r, A) + ^'(r, A) ( 72 ) 

= 2 a;(£+ 1 )+x^ + 2 w 2 (r, A)x(r, A) - x'(u A) ( 73 ) 

= 2 w(£+l) + i? 2 - 2 x'(r,A). ( 74 ) 

Subtracting ( 1711 ) and ( 1771 ) we get 

^'{r,X) = (j)'{r,X) - xir,X) + A, ( 75 ) 
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where A is a constant to be fixed. Integrating we get 

^(r,A) = (f>(r,A)-x(r,A)+Ar + B. (76) 

Both the constants A and B turn out to be zero. This is seen by noting that the large r behaviour 
of the potential (HOI) is correctly reproduced by wi{r,A) term in Wo(r) = wi{r,A) +^(r, A). 
Therefore any addition oi Ar + B with A, B ^ 0 in Wo(r) is not possible, Thus we arrive at the 
result 


Wo{r, A) 


wi{r,A) +<)!<(r, A) - x(r-, A). 

1 ) 


1 {i + 

—ujr - 

2 r 


— Lur 


\9yL^{-y)] 


'dyLt^^\-y) 

[ L^i-y) \ 

+ ujr 



(77) 

.(78) 


In this form the superpotential for the extended radial oscillator agrees with that known in 
literature [3], [TO] . 

The expression for Wo{r,A) can also be obtained by an alternate approach by noting that 
Wo(r, A) is just the superpotential constructed out of the ground state wave function of the 
potential IA(“l(r). The wave functions for V^~\r) are given by (H^ . Taking logarithmic 
derivative of the ground state wave function gives an expression for Wo('r, A) that matches with 
([7R|. The details of this computation are given in the appendix. 

The process of construction of the rational extension of the radial oscillator potential de¬ 
scribed in detail above leads to extended Hamiltonian related to the LI exceptional orthogonal 
Laguerre polynomials. The above process of rational extension of potentials by isospectral shift 
deformation is summarized in the the left half of the flowchart given in figure I. In addition, 
we also give the steps to construct its shape invariant partner in the right half. 


4.4 L2 and L3 rational extensions 

In this section we will briefly summarize our results on rationally extended potentials obtained 
when a different route to the rational extension is taken. 


L2 Category 


First we would like to mention that the entire process leading to the exceptional Laguerre poly¬ 
nomials of Li category can be repeated by exchanging the roles of W 2 {r) and W 3 {r). The two 
extended potentials, obtained by these procedures, correspond to the SUSY partner Hamilto¬ 
nians listed as being related to exceptional Laguerre polynomials of L 2 category by Odake and 
Sasaki [TU] . 

Here, we give the hnal expressions for the L2 exceptional Laguerre polynomials 


{i-m + 1/2) L 


(^ _ ^ + l/2)L-^-3/2(y)i^„+l/2(y) + 


The corresponding weight function is given by 


(79) 


Wm = 


exp(—wr^/d) 


--l-Xjl 


(2/)lu 


(80) 
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Fig.l. 
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L3 Category 


The method of spectral shift deformation when applied to the superpotentials W 2 {t) and w^^r) 
yielded rational extensions of the radial oscillator potential. These potentials were shown to have 
solutions in terms of the LI and L2 Xm exceptional Laguerre polynomials. The above process 
can be repeated with the superpotentials wi (r) and W 4 (r) and we can construct further rational 
extensions of the radial oscillator. From the table it is clear that W 4 {r) = —wi{r, (Z + 1) —?> 1). It 
is seen that the isospectral shift deformation, using these two superpotentials leads to two new 
rational extensions of which one is a singular potential and not interesting in the present context. 
We look at the other potential which is non-singular for certain values of the corresponding 
parameters. In general, it has singularities both in the complex plane and the physical region. 
These singularities are governed by the zeros of the polynomial in the corresponding (j){r). 
The explicit form of the new generalized rational extension and its polynomial solutions are 
presented here for the first time. 

The rational extension obtained by performing the isospectral shift deformation using wi (r) 
and demanding = vj^\r) + R 3 gives 


Vi-Hr) 


\r) + 2 ujr 



^ + i \ dyLZ^i-y) 
2 J L^{-y) 


—uirdy 



dyL^{-y) \ 

L^i-y) ) 


-f 2 



dyL^ii-y) 

L^{-y) 


2 


1 9 


(81) 


where a = —I — 3/2 and again m = 1, 2, 3.... The corresponding eigenfunctions are given by 


where 


Pn,m (^) — 


- _ / ya+^)/^exp(-la;r^) \ 

^n,m{ ) I L“ (ia;r2) ) 


w?’L„“+^(j/)L“ (- 2 /) -h -(m - I - |)L“ ^{-y)L^{y) 
r 2 


are polynomials orthogonal with respect to the weight function 

/+2 exp(-a;rV4) 


W™(r) = 


y^\y—iujr^ 


(82) 

(83) 

(84) 


in the interval 0 < r < 00 . The weight function has singularities governed by the zeros of 
the polynomial L(^(—a;r^/2). From the Kienast-Lawton-Hahn’s theorem on the zeros of the 
Laguerre polynomials [36], [37], the polynomial L(]((—a;r^/2) will not have any zeros on the 
positive real line as the argument of the polynomial is strictly negative. But the polynomial 
will have one negative zero if m is odd and a lies in the range (—m — ^) < a < —m. Thus 
for even m and a not lying in the above given range, the weight function is non-singular in 
the physical range. Thus for all such values of m and a, V^~') (r) has polynomial solutions. 
From the above discussion it is clear that the polynomials Pn,m(j') are orthogonal and form 
a complete set. In addition, the polynomials start with a degree m > 0. Thus we have a 
new set of exceptional polynomials and following the nomenclature in the literature, we name 
them the L3 exceptional polynomials. The discussion of the properties of this new series of the 
polynomials falls outside the domain of this paper and will be taken up else where. 

Thus the isospectral shift deformation using wi{r) leads to a rational extension of the radial 
oscillator, whose solutions are in terms of new set of exceptional polynomials named as the L3 
series. Rational extensions of radial oscillator having solutions other than LI and L2 Laguerre 
EOPs have been discussed in [5] and m- 


5 Trigonometric Poschl-Teller Potential 

The trigonometric Poschl-Teller potential is given by 

V{x) = A{A + l)cosec^xB{B-\-l)sec^ X, A>—^,B>—^. (85) 
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The form of the superpotentials is given by 


w(x,a) = acotx — btanx (86) 

with a collectively denoting the parameters {a, 6}. The possible values of the parameters are 


a = —A — 1, b = B, —B — 1 


and the four solutions of the QHJE correspond to the following choices 


(1) E = -{A + B)^, 

(2) E=-{l + A-B)^, 

(3) E = -{l-A + B)^, 

(4) E=-{2 + A + B)^, 


cr = {A,B} 
cr = {-A-1,B} 
a = {A,-B-\} 
cr = {-A-l,-B -1} 


(87) 


( 88 ) 


The corresponding solutions will be denoted as Wfc, fc = 1,2,3,4, respectively. With the mapping 
r defined as 

t{A,B} = {A+1,B + 1}, 

the shape invariance property is obvious from W 4 ,{x,\) = —wi(x^t{X)) and from wz{x,\) = 
-W2{x,t{\)). 


5.1 Rational Extension of Trigonometric Poschl-Telller Potential 

We note that the first solution in (1551) corresponds to the exact SUSY case and to a normalizable 
ground state solution. Solutions (2) and (3) in (1551) correspond to broken SUSY and will be 
used to construct rational extensions with solutions correspondng to J1 and J2 exceptional 
Jacobi polynomials 0,0, m - 

To keep our intermediate steps general we take 


w{x, a) = acotx — btanx (89) 

and will substitute suitable values for the constants a = {a, b} at the end. The first extension 
process is as follows. 

= w{x)'^ ±w'{x) (90) 

w{x) = w{x,(j) + (j){x), (91) 

where w{x) is determined by demanding 

V^+\x) = V’'+\x) + K, (92) 

where K is a constant. (j>{x) satishes the Riccati equation 

ip^ix) + 2w{x)4>{x) + (j)' {'^) = (93) 

This equation can be linearised by writing 

0(X) = ^ (94) 

u{x) 

and the equation for u(x) turns out to be 

u” + 2w{x, a)u' — Ku = 0 (95) 

u"+ 2{acotx — btanx)u — Ku = t). (96) 


A point canonical transformation y = cos 2x leads to the equation 


,, ,, , , .du 

(1 - 2^ + [b-a + {a + b+ 


Kiu = 0 


(97) 
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where Ki = K/A is a constant. This equation has polynomial solutions with Ki = N{N + 
^ + /j, + 1) with {v, /r} = {a — 1/2, b — 1/2}. With this choice of Ki the polynomial solution 
coincides with the Jacobi polynomial P^'^\y). Therefore, we get 




(dy\ 

rdu{y,X)\ 

\dxJ 

\ dy ) 


— (sin 2x) 


r_^ 

dy 


logP^“ 


- y=cos 2x 


(98) 

(99) 


One can obtain other rationally extended potentials F^“''^(a;) by means of a second extension 
process, as done in the case of radial oscillator. This is done by writing 


w{x) = w{x,a) + x{x) ( 100 ) 

and fixing x(®) by demanding V^~\x) = V^~\x) + which leads to the following equations 
for x{x) 

x{x) = - =-^u" — 2(acotx — bianx)u' — Ku = 0. (101) 

u{x) 

Here again a = {a, b} and can take one of the four possible values listed in (1551) . 

As in the case of radial oscillator, the two superpotentials W 2 (x), W 3 (x) correspond to 
broken SUSY and will give rise to two Shape invarinat nonsingular extended potentials and 
their supersymmetric partners. We obtain two more rational extensions and their partners 
using the other two superpotentials rei(x), W 4 (x) and some of these are likely to be singular. 
Thus starting with the four superpotentials Wk{x), k = 1,2, 3,4, in all, there are eight possible 
extensions obtained by the two processes outlined above and their SI is guaranteed by the 
ansatz that led to equations for the isospectral shift deformation. 

The rational extensions and their partners, obtained from W 2 (x) and wsi^x) will correspond 
to the Hamiltonians related to Jl, and J2 exceptional Jacobi polynomials [3], [IH]. Of the 
remaining four extensions, some will be singular potentials and some others will be regular 
potentials. A detailed study is needed to make a definite conclusive statement whether other 
regular extended potentials are related to the exceptional Jacobi polynomials already found m 
or as in the case of radial oscillator, will lead to new exceptional Jacobi polynomials. Further 
investigation in this direction, along with some other related studies is in progress and will be 
reported elsewhere. 


6 Concluding Remarks 

In this paper we have discussed a procedure to systematically construct rational extensions of 
SIPs. The extended potentials are designed to be shape invariant. Our treatment makes 
it possible to construct an extension which interpolates between various potentials 
related to Xm— exceptional polynomials for different values of m. So for example, 
a solution of (l37ll with arbitrary value of the constant R will give rise to such an 
extended potential, which for special values of R will reduce to potentials related 
to Xm exceptional polynomials for various m values. 

The method can be used to construct rational extenions of all known SIPs [S]. Rational 
extensions of the radial oscillator and the trigonometric DPT have been constructed explicitly. 
In the case of radial oscillator, we obtained the well studied rational extensions with LI and L2 
exceptional Laguerre polynomials as solutions. In addition, we have presented another set of 
infinite generalized rational potentials whose solutions involve a new set of exceptional Laguerre 
polynomials, whose generalized expressions are given for the first time. 

In the case of the DPT potential, we obtain the extended superpotentials from which rational 
extensions of DPT can be constructed. It is clear that this method leads to new rational 
extensions other than those discussed in literature, hence further study is required to see if 
these lead to further new types of exceptional Jacobi polynomials. It would be interesting to 
investigate further generalizations of our method and see where they lead to. 
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APPENDIX 


In this appendix we will give details of and an alternate derivation of (17811 . 

The superpotential and the the Hamiltonian for the radial oscillator are given by 


1 ^ -I- 1 

w{r) = -ujr - 


= 


Z I Z Z I 

p -I- -w r -I- 


£{£ + !) 




' 2 ' 

The energy eigenvalues and the eigenfunctions are 

= 2na;; (y) = exp{-y/2)Li+^/^{y), 


( 102 ) 

(103) 

(104) 


where y = 

The eigenfunctions for as obtained from the above expressions, by replacing i 

£ — 1 are 

= Ny‘^‘^exp{-y/2)L^(y), a =(£-1/2). (105) 

A rational extension was obtained by starting from W 2 = ^ujr + £lr and introducing w(r) = 
W 2 {r) + (j){r). We demanded that (/(r) satisfy U(+)(r) = U(+)(r)-|-i?i and be a rational function 
of r. This led to 

'dyL^i-y)' 


W 2 ,m{r) = W 2 (r) +ujr 


L^{-y) 


(106) 


"01,™ W = (-^+W2.m(r))V'i’^)(r). 

For this purpose first note that of (I105p . which leads to 

d 


= '4’i^\y)x^r 


The eigenfunctions of the extended potential ^ (r) can be computed using the intertwining 
relation 

d ^ 

(108) 

(109) 

( 110 ) 

( 111 ) 
( 112 ) 

(113) 

(114) 


d 
••dr 

dr J \dy 

r_^ 1 , dyL'^iy) 

V2y 


2 K{y) y 

Use of the above result in (11071) gives the following expression for 

^ , 1 dvK^y) 


V'i.mW = V';!:(»')(wr) 


2 ?/ 


2 \{y) 

(a) 


1 £ dyLl^i^-y) 

2 ^ 2y L^{-y) 




LZ,{-y) 9,L“(2/) 


Using the identity 


.■^m(-y) Lniy) 
dxL'^ix) = (x) - L^+^)(x), 


-bl 
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we get 


(115) 


^Ll(r) 


(cjr) 


V L^(-y) 


dyKiy) 

Kiy) 


y—^ujr'^ 




The final expression for the wave function 'tpn^Xir), as obtained by using the equations (I113D 
and (11151) is 


= r^/^exp(-ia;r 2 ) 


L“+1(-2/)L“(2/)-5,L“(2/)L“(-2/) 


L?a{-y) 


(116) 


J y=^u>r^ 


where ol = I — Xjl. This is the result in and the expression for the eigenfunctions can be 
used to have an alternate derivation of (1751) . 


Setting n = 0 in (111511 gives 

/((-) — 


.bi-i ^ f ( y) 




The desired superpotential can be obtained as 

d 


Wo{r) = -^(log^.iW)- 

Substituting the expression in dnzi) in the above equation gives 
Wo{r) = -^(log^i 7 j,(r)) 


n—0 

1 d , —(+), . 

=-log Vo (’’) + 

r dr 


dyL^{-y) dyL^,^^^\-y) 


1 £+1 

= -Lur - 

2 r 


■ ujrdy log 


■^m(-y) 

m-y) 


(a+l) 


i-y) 




which agrees with dZHD. 
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( 120 ) 

( 121 ) 
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